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ABSTRACT 

A  new  method  is  presented  to  explain  how  noise  Is  generated  under  pack  Ice 
by  ridging  of  the  pack  ice.  The  energy  dissipated  during  the  ridging  process 
is  assumed  to  be  the  proper  measure  of  the  noise  source  level.  Noise  source 
levels  generated  by  ridging  are  simulated.  Noise  intensity  at  a  specific  site 
is  calculated  by  summing  signals  from  all  these  sources  after  accounting  for 
propagation  losses.  Calculations  are  made  to  compare  this  simulated  noise 
with  observations  for  an  experiment  conducted  during  the  winter  of  1975-76  in 
the  Beaufort  Sea.  During  a  120-day  period,  461  of  the  intensity  of  the  noise 
signal  is  explained  using  this  process,  and  over  several  20-day  periods,  in 
excess  of  641  is  explained.  In  addition  to  explaining  a  significant  amount  of 
energy  and  ambient  noise,  the  model  is  attractive  on  physical  grounds  and 
properly  explains  lack  of  noise  when  winds  are  high  but  ice  is  strong  enoug’ 

•■.i'  ‘  ••  •’  *.t  ridging. 

INTRODUCTION 

A*  ■  sea  ice  cover  moves,  fractures  and  deforms  it  generates  noise  which, 

In  cot  nation  with  other  background  noise,  interferes  with  sonar,  comraunica- 
tt  <  d  weapons  control  signals  Background  noise  under  sea  ice  has  be,  . 

t  '  undei  a  wide  variety  of  conditions.  Although  various  Invest!  • 
have  i.udied  the  characteristics  ot  this  noise  (so.  Ret  1  r,  ■:>•.),  1 

work  be,  been  done  to  relate  the  observed  noise  to  ape-1  fi,  processes  ’  . 
work  that  has  been  done  includes  attempts  to  relate  ba-»  ground  noise  lev, 
unde-  tea  ice  to  wind  speed,  changes  in  air  temperature.  •  racking  of  th«  1 
and  s  ci.'sse?  and  strains  In  the  Ire  Separate  orrel  -  of  amMeu<  uo 

<1  .  ■  •  •’  sped,  I ,  e  St  re  s.  5  and  ice  strut'-  h.v  .  .  t  '  .  • »  e  ■■  1  o  ’’ 

i  i  hi  i ; :  ,  t  1 1  r  .>  •  ’  i .  i  ■  1 1  "  •  i  •  •  -  • , ,  !  i  •  ■  1  ■  'I 
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noise.  '  *  Better  correlations  have  been  achieved  by  relating  thermal 
cracking  due  to  dally  temperature  drops^*^  and  stress-induced  microcracking. 
These  initial  attempts  to  relate  a  single  noise-making  process  to  ambient 
noise  signals  can  result  in  a  better  understanding  of  and  the  potential  for 
modeling  noise  under  sea  Ice. 

The  single  largest  source  of  noise  in  the  Arctic  in  the  winter  is  thought 
to  be  ridging.  In  thi6  paper,  we  examine  the  relationship  between  observed 
low-frequency  background  noise  and  the  occurrence  of  ridging,  as  simulated  by 
an  ice  model.  Thus,  while  concentrating  on  a  single  noise-producing  process, 
we  can  attempt  to  explain  a  significant  amount  of  the  background  noise. 

The  critical  assumption  made  in  this  work  is  that,  in  ridging,  the  inten¬ 
sity  of  a  noise  signal  at  its  source  (the  ridge)  is  related  to  the  amount  of 
energy  dissipated  in  the  ridging  process.  A  phenomenological  comparison  is 
made  in  which  a  theoretical  estimate  of  noise  due  to  ridging  is  compared  with 
observed  noise  at  a  site  in  the  central  Beaufort  Sea  during  the  winter  of 
1975-76,  The  individual,  small-scale  noise-making  processes  that  occur  in 
ridging  have  not  yet  been  studied.  The  present  phenomenological  study  is  felt 
t'  b  >  necessary  preliminary  effort  in  determining  if  a  more  detailed  study 
von)  '  be  f rui t f u 1 . 

fi  this  study,  a  model  Is  developed  to  estimate  the  not ■;-<  intensity  re 
ceived  at  any  site  from  ridging  at  all  points  on  the  Ice  cover.  This  relation¬ 
ship  depends  on  the  amount  of  energy  dissipated  due  to  noise  production  and 
accounts  for  transmission  losses  during  propagation  The  model  is  applied  t> 
the  Beaufort  Sea  during  the  AIDJEX^  main  experiment  (winter  1975-76)  when 
bacVground  noise  observations  are  available  A  special  method  for  simulating 
•  *  .  developed  to  ns'-  tin  i  ■  'notion  lit  obt  linn.’  froio  r’n  m  mned  can,' 
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and  buoys  deployed  during  that  time.  Data  limitations  restrict  our  study  to 
one-third  octave  bands  centered  around  10  Hz  and  32  Hz.  However,  the  theory 
can  be  extended  to  Include  a  wider  range  of  frequencies.  Comparisons  are  made 
to  estimate  the  fraction  of  total  variance  of  observed  noise  that  can  be 
explained  by  this  model.  Error  estimates  are  provided  to  show  that  the 
unexplained  variance  can  be  attributed  to  data  limitations  and  that  the  basic 
hypothesis  of  the  theory  remains  plausible.  As  a  result,  it  is  suggested  that 
further  work  is  justified. 

I.  NOISE  SIMULATION  MODEL 

The  processes  that  create  noise  can  occur  everywhere  in  the  sea  ice  cover. 
To  describe  the  total  observed  noise  at  a  specific  6ite,  the  sound  waves 
reaching  that  specific  site  from  all  of  the  locations  at  which  noise  is 
generated  must  be  taken  into  account.  The  magnitude  or  intensity  of  the  noise 
at  it  sources  must,  of  course,  be  known,  and  it  must  be  specified  everywhere. 
Practically,  then,  a  model  must  relate  the  noise  source  level  to  some  measur¬ 
able  variable  that  can  be  determined  from  knowledge  of  ice  conditions,  wind, 
current:  temperature  and  other  environmental  parameters.  Finally,  as  the 
arou-'  *  waves  propagate  toward  the  observation  site,  they  are  dispersed  an; 
at  i«i’  <?ed  geometrically  by  reflection  off  the  sea  floor  and  the  underside  o* 
the  t  e  cover  and  by  sound  speed  gradients.  The  approach  presented  her. 
attempts  to  account  for  all  of  these  factors,  thereby  providing  a  mathematical 
model  formulated  on  basic  physical  principles  of  acoustics. 

Many  different  processes  occur  that  generate  noise.  The  types,  of  process.- 
and  cln  relative  importance  of  their  contributions  to  the  total  observed  ooi-;< 
by  bo  -  h  season  and  location.  Tb  I  w-i’  .■'lii.-i'n'  rnffr-  win!  t-r.  i  o  i 
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the  central  Beaufort  Sea.  In  these  heavy  ice  conditions,  ridging  is  generally 
believed  to  be  the  single  largest  noise  source.  But  even  though  ridging  i6 
often  thought  of  as  one  large-scale  geophysical  process,  it  is  in  fact  a  com¬ 
bination  of  many  small-scale  processes.  During  ridge  building,  sea  ice  i6 
fractured  in  bending,  ice  floes  slide  against  one  another  to  be  crushed  and 
sheared,  and  the  resultant  blocks  are  piled  up  and  down  against  gravitational 
attraction. 

A.  Geophysical-scale  model 

Sea  ice  dynamics  understanding  has  increased  immeasurably  in  recent  years, 
and,  fortunately,  the  capability  to  simulate  sea  ice  behavior  has  kept  pace. 

g 

Extensions  to  the  A1DJEX  model,  on  which  the  present  work  i6  based,  have 
been  reported  in  Ref.  9.  Of  these  extensions,  the  relationship  developed  to 
estimate  the  mechanical  energy  budget  of  the  ice  cover^  is  the  most  critical 
for  this  work.  Here  it  is  hypothesized  that  energy  dissipated  by  large-scale 
failure  processes  such  as  ridging  is  the  proper  measure  of  energy  available 
for  conversion  into  acoustic  signals.  Thus,  the  capability  to  estimate  these 
dissipation  terms  Is  necessary  to  the  success  of  the  present  work,  Mori 
detailed  results  on  these  modeling  capabilities  have  been  presented  els 
vh  m e  and  estimates  have  been  made  of  the  Beiufort  Sea  energy  budget 
during  AIDJEX . 1 ' 

The  small-scale  mechanisms  in  ridging  dissipate  energy  that  provide 
energy  sources  for  acoustic  signals.  If  these  small-scale  mechanise,  energ. 
•inks  are  to  be  useful  for  explaining  ambient  noise,  there  must  be  a  way  i 
measure  the-  at  different  times  and  locations  Fortunately,  the  AID-Jt-X  m ■  ■ 
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does  provide  a  vay  to  calculate  these  energy  sinks  in  terms  of  the  geophysical- 
scale  behavior.  Specifically,  in  a  constitutive  law,  it  is  assumed  that  the 
energy  dissipated  by  the  geophysical-scale  stress  during  deformations  is  equal 
to  the  energy  dissipated  by  these  small-scale  mechanisms.  In  the  model, 
expressions  are  available  to  estimate  these  energy  dissipation  levels  as  a 
function  of  ice  conditions  at  different  times  and  locations. 

The  large-scale  constitutive  behavior  is  represented  by  an  elastic-plastic 
mathematical  model.  A&  in  other  plasticity  models,  the  stress  is  limited  to 
lie  within  a  yield  surface.  However,  a  stress  resultant  o  is  used.  This 
stress  resultant  is  the  Cauchy  stress  integrated  through  the  ice  thickness. 

The  theory  is  two-dimensional. 

Ar  tsotropic,  diamond-shaped  yield  surface  is  used  in  our  computations. 

Ih  yield  criterion  is  given  by 


<J> (a x ,0 x j.)  <  0 


(1) 


where  |  tr  g  and  ■  -j  tr  with  o*  ■  O  -  Oj  1.  The  invariants 

Oj  and  Ojj  are  half  the  sum  and  difference  of  the  principal  stress  value, 
respectively,  and  p*  is  the  ice  strength  in  isotropic  compression.  The  shear 
stress  Invariant  0^  is  nonnegative  when  defined  in  this  way.  Plastic 
deformations  are  assumed  to  satisfy  an  associated  flow  rule  Thus,  the 


pla  t.  ( .  stretching  is  orthogonal  to  the  yield  surface  and  may  be  expressed  as 


D 

-r 


(2) 


where  A  is  a  nonnegative  scalar, 

'  linear  elastic  law  is  assumed  when  stiess  is  within  the  yield  surfact 

£  *  (Mj  Hj)l  ir  f  i  ? 
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where  e  la  elastic  strain,  and  and  M2  are  the  bulk  and  shear  moduli,  respec¬ 
tively.  The  elastic  strain  Is  related  to  stretching,  D,  spin,  W,  and 

13 

plastic  stretching  by  the  relationship 

e-We  +  eW-D-D  .  (4) 

~P 

Ice  strength  Is  determined  as  a  function  of  the  ice  thickness  distribu¬ 
tion.  Changes  in  the  strength  occur  as  deformations  change  the  thickness 
distribution,  which  Is  governed  by  the  relationship 

G  +  ft  ||  -  ~  GV-v  (5) 

where  G(h,t)  is  the  fraction  of  area  covered  by  ice  thinner  than  h  at  time  t, 

and  4*  Is  the  rate  of  ice  production  of  ice  thinner  than  h  due  to  mechanical 
14 

redistribution,  Assuming  4>  Is  linear  in  the  rate  of  plastic  stretching,  then 

-  D  [a  (0)  +  a  (6)W  )  .  (6) 

p  o  r  r 

Here,  Dp  -  f(Dp)  i;2  +(Dp)lI2]1/2  where  (Dp)  l  and  (Dp)  H  are  the  sum  and  difference  of 

principal  stretching  values.  The  angle  8  *  arctan  ^Dp  )ll/(Dp)l]  16  the  ratio  of 

shearing  to  dilating.  The  coefficient  aQ  is  the  fraction  of  open  water  created 

by  a  unit  of  deformation.  Conservation  of  mass  requires  a  and  a  to  be  related 

or 

by  •  a  *  cos  P.  The  variable  W  describes  the  conversion  of  thin  ice  int> 
or  r 

thick.,  Ire  by  ridging  and  is  3  complicated  functional  of  the  thickness  di<='ri- 
buc 1 o  ’ 

Th  rate  at  which  energy  is  dissipated  by  large-scale  stresses  during 
def  ormation  is  determined  as  the  product  of  stress  and  permanent  stretching 
This  1  k  the  dissipative  stress  power,  p ^  *■  tt  o  D  This  large-scale  meas  >r.  of 
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nergy  dissipation  is  equated  to  the  sun  of  identifiable  snail-scale  energy 


inks.  Thus, 


rhere  is  the  set  of  small-scale  energy  sinks  that  are  known  to  dissipate 
nechanlcal  energy.  This  relationship  provides  an  estimate  of  the  large-scale 
Ice  strength  and  helps  relate  the  redistribution  function  to  the  constitutive 
law.  For  this  study  of  noise  sources,  the  relationship  also  provides  a  means 
of  estimating  the  magnitude  of  the  noise  source  level.  We  assume  that  each 
small-scale  energy  sink  is  a  source  of  noise  and  that  the  energy  dissipated  by 
each  sink  can  be  determined  from  the  large-scale  energy  dissipated. 

Gravitational  potential  energy  changes  and  dissipation  by  friction  as  ice 
hloiks  slide  through  the  keel  and  sail  are  considered  to  determine  the  force 
neerft  to  build  a  ridge. This  concept  has  been  generalized  by  allowing 
both  compression  and  shear. ^  The  change  in  gravitational  potential  energy, 


,  is  ice  is  piled  up  and  down  into  the  ridge  sail  and  keel  is  given  by 


q  -  c 
P  P 


OC 

/ 


k.  «  i 

P  2 


.he  gravitational  acceleration,  and 


D  0  (Cw  '  f  ' 
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In  concept,  more  buoys  would  allow  better  resolution  and,  thus,  would  reduce 
this  potential  error.  However,  the  sea  ice  constitutive  law  assumes  that  many 
ridges  are  forming  within  a  triangle  to  ensure  isotropy  and  homogeneity,  bo 
10  km  is  roughly  the  lower  limit  for  achieving  better  resolution  with  a  more 
dense  buoy  array.  In  addition,  as  buoys  approach  one  another,  the  uncertainty 
in  stretching  from  the  velocity  inhomogeneity  variation  becomes  so  large  that 
it  dominates  errors.  Therefore,  triangles  roughly  on  the  order  of  50  km 
provide  data  that  are  as  accurate  and  consistent  as  is  reasonable  for  using 
this  simulation  technique. 

Another  potential  source  of  error  is  due  to  the  assumed  thickness  distri¬ 
bution  used  in  the  simulations.  The  initial  thickness  distribution  is 
reasonably  accurate  for  Julian  day  465,  It  is  used,  however,  to  begin  the 
simulations  for  each  20-day  block.  This  choice  affects  the  strength  level, 
which,  if  constant,  would  not  affect  the  comparison,  and  it  affects  variations 
in  strength  due  to  redistribution.  The  effects  of  this  assumption  are  not 
known.  Our  inability  to  measure  the  thickness  distribution  of  sea  ice  is  a 
problem  with  all  simulations  of  sea  ice  behavior. 

!  CONCLUSIONS 

In  this  paper  we  have  introduced  a  new  concept  to  explain  background  noise 
under  Arctic  sea  ice.  The  concept  is  based  on  the  hypothesis  that  energy 
dissipated  by  the  ice  in  ridge  building  is  related  to  the  intensity  of  the 
source  level  of  background  noise.  A  mathematical  model  of  sea  ice  behavior  is 
capable  of  simulating  the  energy  dissipated  by  stress  during  deformation 
everywhere  In  the  i  c  covt»  r  The  spatial  distribution  of  noise  sou  eves  idj  * 

•*!’■>*  *  h*  s  1  n  il »  c(*  a  t  h  n  v  tiwr  It  f.  r«ins«*.f  ss  i -v  loss'"'  ar.  *»:  i">T‘ 
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as  a  typical  value  for  the  time  period  used  in  our  analysis.  This  uncertainty 

in  the  variation  in  velocity  gives  rise  to  uncertainties  in  stretching  that  are 

0.7  x  10  ^  s  *  at  the  90%  confidence  level  during  spring.  These  errors  assume 

a  gauge  length  of  100  km.  Smaller  scale  averages  generate  proportionately 

larger  uncertainties  in  stretching.  Since  triangles  are  typically  SO  km  in  6ize, 

uncertainties  in  stretching  are  on  the  order  of  1.5  x  10  ^  s  *. 

It  is  difficult,  if  not  impossible,  to  estimate  uncertainty  in  the  power 

dissipated  because  the  extremely  large  uncertainty  in  stretching  makes  uncer- 

tainty  in  stress  unrealistically  large.  For  a  strength  of  5  x  10  N/m,  the 

-3  2 

stretching  uncertainty  gives  rise  to  an  uncertainty  of  7  x  10  W/m  .  If  this 

is  input  as  a  constant  error  over  a  circular  domain  with  a  radius  of  500  km  and 

the  transmission  loss  function,  Eq.  (22),  is  assumed,  then  integration  gives  an 

9 

error  ir  simulated  noise  of  about  1  x  10  W.  This  theoretical  maximum  error 
is  js  targe  as  maximum  simulated  noise  values  (Figures  3  and  4).  Another 
estimate  of  velocity  Inhomogeneity  errors  is  obtained  by  looking  at  quiet  times 
when  all  simulated  noise  may  be  assumed  to  be  due  to  velocity  inhomogeneities. 
These  values  were  discussed  earlier.  Fortunately,  the  actual  values  are  about 
25%  tf  the  simulated  maxima,  far  lower  than  the  theoretical  maxima. 

Uncertainties  in  the  transmission  loss  law  also  limit  the  aevuraev  of  tills 
sim>1aL,.'n  technique.  First,  the  transmission  loss  law  is  a  emj  i  :  i  al  fi>  to 
data,  and  some  unexplained  variance  must  by  expected.  Second,  the  cransni  ss !  •'« 
los*  law  decays  exponentially  from  the  source.  Thus,  any  sources  very  n<  n  the 
observation  location  will  be  very  sensitive  to  the  location  of  the  soum  Bat 
a  sourr-  can  be  located  to  only  be  within  a  specific  triangle,  rather  than  a’ 
an  »*.»■  site  within  the  triangle.  This  Is  i  limitation  of  the  basic  app;  >■<. 
i  'n  d  1  s’  i  •  •  ■  be f  v»  c  ,  bo-  (  i  t  1 1  i  i  <-  •  m  t  r  - . !  ; ,  '  h  r  t  '  nc  ‘ 
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Scatter  plots  of  observed  signals  and  simulated  noise  signals  from  qr  are 
given  in  Figures  5  and  6  for  10  and  32  Hz,  respectively.  These  plots  show  that 
the  observed  and  simulated  signals  tend  to  rise  and  fall  together.  Many  signals 
occur  at  low  noise  levels,  which  causes  a  problem  in  the  plots  because  the  data 
points  tend  to  black  out  the  area.  It  would  have  been  simple  to  plot  data  at 
logarithmic  scales  to  spread  out  the  points  at  low  noise  levels,  but  this  would 
have  emphasized  data  that  are  subject  to  relatively  large  errors.  Several 
individual  events  appear  in  the  scatter  plots.  These  have  been  Bhown  by  con¬ 
necting  the  dots  to  make  a  "spaghetti*'  plot.  The  path  taken  by  the  larger 
events  clearly  shows  the  correlation.  Some  difference  appears,  though,  86 
individual  events  are  related  by  different  regression  equations  or  slopes. 

This  is  probably  caused  by  our  lack  of  resolution  (triangle  size),  by  inaccura¬ 
cies  in  the  transmission  loss  law,  or  by  our  uncertainty  in  ice  conditions  and 
therefore  Ice  strength  during  the  experiment . 

The  data  are  obviously  shifted  to  the  right  and  do  not  pass  through  the 
origin  This  is  caused  by  excessive  energy  being  dissipated  as  individual 
buoy;  respond  to  small-scale  velocity  inhomogeneities  that  are  not  modeled  by 
ihe  3 ci  ice  constitutive  law.  If  the  ridging  process  were  modeled  perfectly, 
one  wr.ild  expect  the  regression  line  to  Intersect  the  ordinate  and  not  th 
absciss.*  ,  because  other  noi se -mnki ng  processes  woild  generate  noise  thy  h.i 
n >t  been  modeled  here. 

Although  errors  can  occur  because  spatial  variations  In  the  velocity  field 
are  noi  smooth,  the  large-scale  average  is  a  reasonable  estimate  of  the  deforma¬ 
tion.  occurring  within  the  triangle.  The  inhomogeneity  variation  has  beer 

22 

S- *  fm.it.  ■■{  to  be  on  the  order  of  0 ,  '*  cm/s  during  spring  1975.  We  use  th! 
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the  comparison  varies  throughout  the  4-month  period.  In  addition,  three  energy 
dissipation  variables  are  compared  with  observed  noise  to  determine  which  one 
explains  the  most  variance  In  observed  noise.  For  the  entire  period,  the  energy 
dissipated  by  ridging,  qr,  has  correlation  coefficients,  r,  of  0.67  and  0.68  at 
10  and  32  Hz,  respectively.  The  energy  dissipated  by  shearing,  qg,  explains 
much  less  of  the  variance,  with  correlation  coefficients  of  0.36  and  0.37  at  10 
and  32  Hz,  respectively.  Both  coefficients  vary  greatly  from  one  20-day  block 
to  another.  For  completeness,  p^  Is  also  used  as  a  measure  of  energy  dissi¬ 
pation  to  simulate  noise.  It  does  not  do  as  well  as  q^,  except  at  times  when  q 
Is  better,  because  p.  »  q  +  q  and  both  q  and  q  are  correlated  somewhat. 
During  Julian  days  361  through  380,  ridging  explains  very  little  of  the  noise, 
whereas  shearing  explains  a  larger  amount.  This  is  a  rather  quiet  time,  but 
one  during  which  energy  is  dissipated  in  numerous  high-frequency  spikes.  For 
Julian  days  381  through  460,  the  correlation  coefficient  is  consistently  higher 
for  noise  simulated  from  q^.  However,  during  the  period  between  Julian  days 
461  and  480,  q  is  again  greater  than  q  ,  but  for  both  the  correlation 
coefficients  are  high.  The  threshold  in  simulated  noise  reduces  the  correlation 
when  noise  levels  are  low.  Lag  correlations  were  also  determined.  Maximum 
correlations  appear  at  zero  lag.  Correlations  are  about  as  high  at  lag  times 
(<f  1  h.  but  not  substantially  higher.  This  is  likely  caused  by  the  original 
filtering  of  motion  and  noise  data  to  remove  high-frequency  inpwf  To  learn 
further  how  well  simulated  noise  can  explain  observed  noise,  a  regression 
analysis  was  performed.  Multiple  regressions  of  observed  noise  onto  eneig 
dissipated  by  ridging  and  shearing  do  not  explain  a  more  appreciable  annum  of 


?!<♦  variance. 
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locations  because  they  are  caused  by  atmospheric  pressure  systems  moving  across 
the  Beaufort  Sea.  One  can  expect,  then,  that  many  variables  vlll  rise  and  fall 
at  roughly  the  same  time. 

The  simulated  noise  Is  nearly  always  nonzero,  while  the  observed  noise  drops 
to  nearly  zero  at  times.  This  threshold  for  the  simulated  noise  is  attributed 
to  velocity  lnhomogeneitles.  These  uncorrelated,  random  motions  of  Individual 
Ice  floes  dissipate  energy  In  the  model  because  of  the  dissipative  nature  of 
plasticity.  An  estimate  of  the  size  of  this  contribution  may  be  made  during 
two  quiet  periods,  Julian  days  411  through  421  and  451  through  461.  The 
10-  and  32-Hz  simulated  noise  levels  from  have  values  on  the  order  of 

o  11 

0.2  x  10  and  1.0  x  10  mW/m  ,  respectively.  The  threshold  levels  for 

noise  from  q  are  about  half  these  values, 
s 

On  a  few  occasions,  either  the  simulated  or  the  observed  noise  curve 
indicates  energy  input  but  the  other  does  not.  Noise  observed  in  the  absence 
of  a  simulated  noise  signal  can  occur  if  it  is  generated  by  a  process  other 
than  ridging.  On  the  other  hand,  simulated  noise  can  occur  in  the  absence  of 
observed  noise  if  the  model  dissipates  energy  at  a  location  too  close  to  the 
sensot .  The  simulated  noise  is  very  sensitive  to  the  location  of  energv 
dissipation  when  near  the  sensor,  i.e.,  within  abour  100  km.  Since  thf 
location  of  the  dissipated  energy  can  only  be  resolved  to  within  abou*  50  k " , 
thin  limits  the  accuracy  to  which  the  nearby  noise  can  be  simulated. 

8  Coi relation 

A  quantitative  measure  of  comparison  Is  obtained  by  determining  the  cor 
relation  between  observed  and  simulated  noise.  Correlation  coefficients  ar. 

d  in  Tahir  II  Tin-  t  Ini'  si-i  i  •  Is  bred  i  n  into  ?•"  <1  i  bio-.-  t  •  1  •  •  it 
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data  lor  the  10-  and  32-Hz  noise  appear.  The  observed  noise  and  the  ice 
velocity  data  have  been  filtered  to  remove  variations  vith  a  period  shorter 
than  one  day.  All  data  are  sampled  at  synoptic  time  Intervals.  Simulated 
noise  has  been  presented  using  both  the  energy  dissipated  in  ridging,  qf,  and 
the  energy  dissipated  In  shearing,  qg  as  energy  dissipation  sinks. 

These  time  histories  show,  in  a  qualitative  way,  the  ability  of  the  model 
to  simulate  background  noise.  Intensities  of  the  acoustic  signals  are  pre¬ 
sented  on  a  linear  scale.  Observed  noise  is  measured  in  units  of  square  micro- 
pascals  per  hertz,  (pPa)  /Hz.  This  is  the  power  of  the  acoustic  wave  per  unit 
area.  The  values  for  the  simulated  noise  are  given  in  milliwatts  per  square 
meter,  mV/m  .  This  scale  is  arbitrary  because  the  signature  function  S(f)  is 
not  included  in  the  calculation.  Since  it  is  a  constant  for  each  frequency  band, 
a  comparison  may  be  made  in  the  form  shown,  A  linear  regression  Is  used  to 
provide  an  estimate  of  the  value  of  the  signature  function  at  each  of  the  two 
freq  'ancles . 

A  visual  comparison  of  the  observed  and  simulated  noise  from  qf  is  good. 

However,  the  comparison  between  observed  and  simulated  noise  from  qg  Is  not 

as  good.  Both  observed  and  simulated  time  histories  are  similar  for  the  10  and 

32  '1/  signals.  The  variations  generally  occur  at  the  same  times  and  have 

sla*.  '.  «<  durations.  Furthermore,  the  sizes  of  the  larger  events  are  comparable. 

Thl  omparison  is  especially  good  during  the  most  intense  events  (for  example, 

between  days  392  and  411  and  days  438  and  443).  Activity  appears  to  rise  ar 

roughly  5-day  intervals  as  storms  pass  through  the  Beaufort  Sea.  This  relates 

1  2 

to  the  energy  input  by  the  wind  field  Previous  estimates  of  the  energy 
I  op  by  winds  and  currents  during  the  period  under  consideration  her.  provid* 

.• :  n.  «  result'.  Thor,  f »  i  strung  rot  r.'l  at  !  on  betw.or  win!  a'  diri •  r  •  •  n 
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of  sinulated  noise  require  that  the  integrand  in  (25)  be  evaluated.  We  have 
assumed  that  the  power  dissipated  is  homogeneous  within  each  triangle  so  that  it 
may  be  removed  from  under  the  Integral.  The  transmission  loss  remains  to  be 
Integrated.  A  numerical  integration  scheme  is  chosen  using  values  of  the  inte¬ 
grand  at  the  midpoints  of  the  sides  (denoted  by  subscripts  a,b  and  c)  so  that 

f  T(|x-*|,s,f)  da  -  j  [l(lRal,s,f)  +  TdRj.s.f)  +  T(|Rcl,s,f)J  .  (26) 


This  is  a  simple  formula  that  is  accurate  to  the  second  order. 

A  higher  order  integration  formula  could  have  been  used  to  Integrate  the 
transmission  loss  contribution  within  each  triangle.  Such  accuracy,  however, 
would  ignore  the  fact  that  all  variations  of  q  have  been  neglected  and,  there¬ 
fore,  would  not  provide  increased  accuracy  in  the  simulated  noise.  Also,  the 
simplest  formula  expressed  in  terms  of  values  at  the  corners  is  troublesome 
when  considering  those  triangles  at  whose  corners  the  noise  is  being  simulated 
because  the  transmission  loss  is  not  appropriate  for  propagation  ranges  less 
than  i  7  km  (20  nautical  miles). 


I  i  -  RESULTS 

Background  noise  simulated  using  the  model  is  compared  with  observed  noise 
for  one  site  in  the  Beaufort  Sea  during  a  4-month  period  in  the  winter  ot 
1975-76.  A  time  history  of  data  is  presented.  Statistical  correlations  anti 
line. a i  regressions  are  used  to  quantify  the  comparison 


A  \ \  >  hist  or  y 

'lirert  oor,|  a  r  i  son  between  observe!  and  Simula’ <!  *  i  historic 

«ph<»n<  si'  'itin’v  1'  (Figur*  2 )  is  pre-in*.  1  ,  -  Fi.  ■  •  •  1  i 
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When  the  array  was  deployed,  the  Beaufort  Sea  region  was  subdivided  into 
38  triangles  as  shown  in  Figure  2  (solid  lines).  After  Julian  day  421,  buoy 
number  16  moved  rapidly  Into  the  Chukchi  Sea,  and  these  triangles  became  too 
distorted  for  meaningful  simulations;  the  five  triangles  surrounding  the  buoy 
collapsed  into  three  triangles  (dashed  lines). 

The  Integral  In  Eq.  (24)  sums  noise  signals  from  surrounding  locations.  It 
may  by  evaluated  as  a  sum  of  integrals  over  the  set  of  triangles.  The  advantage 
of  introducing  the  triangles  is  that  the  energy  sink  q  may  be  estimated  simply 
in  each  smaller  triangular  region.  Since  q  is  calculated  as  the  product  of 
stress  and  stretching,  both  variables  must  be  calculated  in  each  triangle. 
Average  stretching  is  estimated  directly  from  velocities  of  buoys  located  at 
each  corner.  Stress  is  estimated  by  solving  the  constitutive  law  by  inte¬ 
grating  the  stretching  history.  The  product  then  gives  the  desired  energy 
dissipation . 

The  motions  of  three  buoys  are  used  to  determine  the  velocity  gradient, 

L  ■  Vv.  The  average  may  be  found  by  using  the  Green-Gauss  theorem  to  express 
the  integral  as  a  line  integral  and  then  assuming  linear  variation  of  velocity 
alont  the  boundary.  The  relationship  becomes 

A  L  *  /  v  ®  n  d£  (25) 

*S  ~ 

where  L  is  the  gradient  of  the  average  ice  velocity,  v,  over  the  region  R 
bounded  by  the  curve  S. 

fhe  stress  history  is  obtained  from  the  elastic  plastic  constitutive  law 

while  assuming  homogeneous  stretching  within  each  triangle,  Equations  (11 

20 

through  (5)  are  integrated  numerically  using  a  difference  scheme  as  part 
■  '  mirac-rical  scheme  to  solve  the  sen  («•■  dynamics  moth  '  Actual  os  t  :  tr.a :  i" 
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Therefore,  If  we  are  able  to  estimate  the  signature  function,  S(f),  for  a  selec¬ 
ted  noise-making  process,  then  the  noise  that  will  be  observed  at  a  site  x 
from  the  process  will  be  given  by  (24)  in  terms  of  the  energy  dissipated  at  all 
locations  and  times.  We  have  included  the  time  variable  t  to  show  explicitly 
how  temporal  energy  variations  In  energy  dissipation  cause  time  variations  in 
noise  generation. 

11.  SIMULATION  TECHNIQUE 

During  the  AIDJEX  main  experiment  (the  winter  of  1975-76),  an  array  of  buoys 

18  19 

and  omnidirectional  hydrophones  was  deployed  in  the  Beaufort  Sea.  ’  The 
hydrophones  were  suspended  at  a  depth  of  100  ft  below  the  top  of  the  ice.  The 
pressure  measurements  were  converted  to  acoustic  intensities  using  a  power 
meter  They  were  band-pass  filtered  into  one-third  octave  bands  centered  on 
3.2,  10,  32  and  1000  Hz  (in  one  case,  100  Hz  replaced  32  Hz).  Signals  were 
averaged  over  45  s  and  sampled  at  3-hr  synoptic  time  intervals.  The  3.2-Hz 
data  were  contaminated  by  noise  from  cable  strum.  The  other  data  are  felt  to 
be  adequate,  although  sampling  frequencies  are  too  low. 

Figure  2  shows  the  spatial  distribution  of  the  AIDJEX  sensors  on  Julia  day 
38j  (January  16,  1976)  .  The  data  are  adequate  to  estimate  dissipation  b 
ridgl  i  ind  to  compare  simulated  noise  with  the  levels  observed  at  one 
hydrophone  site,  site  number  11.  To  account  for  all  noise  received  at  a  site, 
sound  waves  propagating  to  it  from  all  directions  must  be  included  in  th, 
model  This  requires  that  ice  motion  data  be  available  all  around  the  sit  . 
whic'  limited  study  to  site  number  11,  the  only  hydrophone  site  in  the  cen:. 

■■■>  I  'i'loy  array 
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As  the  acoustic  signal  propagates  from  the  source,  it  undergoes  geometrical 
dispersion  and  dispersion  due  to  reflections  off  the  sea  floor  and  the  underside 
of  the  ice  cover  and  to  gradients  in  the  water  column.  An  empirical  model  has 
been  presented  for  long-range,  low-frequency  (i.e.,  below  100  Rz),  deep-water 
(greater  than  1000  m)  transmission  loss,  TL,  in  the  Arctic. ^  The  equation  16 

TL  -  a  +  10  log  R  +  bf  +  csR  +  dfsR  (22) 

where 

R  is  the  range  in  kilometers, 

s  is  the  standard  deviation  of  ice  depth  in  meters,  and 
f  is  the  frequency  in  hertz. 

The  constants  a,  b,  c,  and  d  are  given  In  Table  I  as  a  function  of  path  length. 
The  loss  function  T  16  related  to  transmission  loss  TL  by 

T  -  10-TL/1°  .  (23) 

Figure  1  presents  the  loss  function  for  waves  with  frequencies  of  3.2,  10,  32, 
and  100  Hz.  The  transmission  loss  law  given  in  Eq .  (22)  is  evaluated  for 
frequencies  from  3.2  to  100  Hz.  Each  curve  has  been  normalized  to  unity  at  a 
rang.:  of  30  km.  Within  this  range  the  relative  intensity  is  assumed  constant. 
Tbi  .  assumption  is  compatible  with  the  restriction  that  the  law  is  valid  onl_. 
fo'  propagation  paths  in  excess  ot  37  km  (20  nautical  miles).  It  will  be  seen 
that  thi6  assumption  is  also  comparable  with  the  resolution  of  ice  motion  data 
used  in  this  study  to  simulate  noise. 

11  the  assumed  relationship  between  the  source  levt?l  noise  and  dissipate  l 
energy  Ik  substituted,  we  find  that 

N(x,f,r)  -  S(f)  J  q(v,t)  T(  !x-yl  ,r.,f)  da  (24) 

A 
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It  Is  worth  noting  that  q  and  therefore  Ng  depend  explicitly  on  the  horizontal 
position  3.  We  do  not,  however,  attempt  to  describe  the  dependence  on  depth, 
being  content  at  this  time  with  a  two-dimensional  theory,  because  data  are 
available  at  only  one  depth  in  our  data  set.  Other  variables  such  as  ice 
conditions  affect  only  q  and  do  not  appear  in  any  other  ways  in  the  nol6e 
source  definition. 

C.  Simulated  noise 

At  each  location  at  which  noise  is  generated  it  is  necessary  to  sum  source 
level  noise  contributions  from  all  noise  generation  mechanisms.  If  it  is 
assumed  that  there  are  enough  sources  that  the  acoustic  signals  are  uncor¬ 
related,  then  waves  must  cause  interference,  but  Intensities  may  still  be 
added  together  linearly.  This  is  a  consequence  of  the  fact  that  intensity  is 
the  iquare  of  the  pressure  in  an  acoustic  signal  and  that  uncorrelated 
acoustic  waves  are  Included  in  a  mean  square  sense. 

The  intensity  of  noise  at  any  site  Is  obtained  by  adding  together  the  con¬ 
tributions  from  all  surrounding  sources.  This  Is  accomplished  by  integrating 
over  the  region  of  interest.  A,  with  the  Integrand  being  the  propagated  noise 
Intensity  ppr  unit  area.  Thus,  we  write 

N  (x,t)  -  J  N  (£,f)  T(  Jx-jrl  ,s,f )  rii  (21) 

A 

where  x  is  the  site  at  which  the  noise  is  to  be  predicted  or  observed,  y  is  the 
integrand  in  the  generic  location  of  each  point  in  the  region,  |x-y!  is  the 
range  of  the  propagation  path  and  T  U  a  loss  function  (the  relative  intensity 
of  th*1  acoustic  wave  after  propagating  from  y  to  x) 
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B.  Noise  source  Intensity 

Consider  the  key  assumption  that  small-scale  energy  sinks  associated  with 
large-scale  ridging  deformations  are  related  to  the  source  level  Intensity  of 
acoustic  signals  generated  by  these  mechanisms.  This  assumption  Is  attractive 
for  two  reasons.  First,  both  quantities  are  measures  of  energy  and  thus  this 
relationship  Is  more  physically  meaningful  than  a  relationship  between  energy 
and  another  variable  such  as  strain,  stress,  ice  velocity,  wind  6peed  or  air 
stress.  Second,  several  limiting  cases  are  identified  where  the  energy 
relationship  is  sensible  and  other  relationships  are  not.  These  include 
situations  in  which 

(a)  Internal  ice  stresses  are  large,  but  not  large  enough  to  cause 

plastic  deformation;  in  this  case,  no  energy  is  dissipated  by  ridging 
and  no  noise  is  expected. 

*b)  Strains  are  large  but  stresses  are  small  because  of  a  large  fraction 
of  open  water;  in  this  case,  only  a  little  energy  Is  dissipated  and 
little  noise  is  expected. 

The  simplest  assumption  that  can  be  made  Is  that  the  intensity  of  each 
noise  source  is  linear  in  the  rate  of  dissipation  of  energy  by  that  mechanise. 

Tiechanistn  is  expected  to  produce  noise  In  a  different  range  of  freque 

cte  Therefore,  a  scalar  signature  function  relating  the  power  dissipau-d, 

q,  t.o  the  noise  source  level  is  introduced.  This  signature  function,  S(f), 

•us*  depend  on  frequency.  Each  mechanism  will  have  its  own  signature  or 

dependence  on  frequency.  If  the  amplitude  of  the  noise  source  intensity 

generated  by  the  mechanism  is  N  (f),  then  the  noise  source  level  is  giver. 

s 


for  a  typical  process  as 


N  (  x  ,  f  '  i !  *' )  q  (  x  ' 

C  — 
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In  addition  to  defining  strength,  relating  dissipative  stress  power  to  energy 
dissipated  by  small-scale  sinks  also  Introduces  a  constraint  between  the  yield 
surface  shape  and  the  coefficient  tnat  controls  the  amount  of  open  water 
formed  lr.  compression  and  shearing. 

The  energy  dissipated  in  shearing,  qfi,  is  assumed  to  depend  linearly  on 

9 

the  rate  of  permanent  stretching  and  on  ice  strength  : 

qs  "  Dp  °s  p*  <19) 

where  the  coefficient  o 6(0)  is  determined  from  the  coefficient  ar  and  the  yield 
surface. 

This  sink  has  not  been  analyzed  by  considering  a  small-scale  mechanism  but 
has  been  introduced  as  a  residual  in  the  large-scale  model  energy  budget.  The 
shear  energy  sink  describes  energy  dissipated  without  redistribution  of  ice. 

It  is  possible  to  think  of  this  term  as  representative  of  energy  dissipated  in 
building  shear  ridges.  Since  shear  ridges  are  narrow  features,  they  do  not 
affect  the  thickness  distribution  significantly.  However,  the  size  of  the 
shear  energy  term  is  determined  without  consideration  of  any  specific 
mechanism.  When  the  shear  sink  is  included,  large-scale  strength  still 
satisfies  (19),  but  there  is  no  longer  a  constraint  between  the  yield  surface 
shap?  and  the  formation  of  open  water.  Instead,  both  variables  are  specified 
independently  using  observations  and  results  of  simulations  of  ice  behavior. 

Integration  of  the  thickness  distribution  equation  requires  that  an  initial 
condition  be  described.  This  quantity  can  have  a  strong  impact  on  the  nois 
simulation  because  It  affects  stress.  We  have  chosen  an  Initial  thicknes 
d ! ii rlbut ion  observed  by  investigators  on  the  USS  Gurnard  during  a  cruis' 
i»r-!>  t!p  Beaufort  Se  t  Ire  during  Api  i  1  19?c. 
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where  p  Is  the  density  of  the  Ice  and  that  of  the  water.  If  the  redis¬ 
tribution  function  (6)  is  substituted,  then 


q  *Da  p* 
P  P  r  *P 


where 


p"  *  c 
P  P 


w 

/,  3W 

"V 


dh  . 


(11) 


(12) 


The  power  dissipated  by  friction,  q^ ,  has  the  same  functional  fora  as  gravita¬ 
tional  potential  energy  changes  and  may  be  written  as 


where 


qf  -  Dp  cr  p* 


Cf  * 


.  o  f  h2  a(h) 

cf  /  r-(l7k)  dh 

*,0 

P'(P*  ~  P>»  [P(k-l?|  ' 
2  tan  $  *  [  Pw  J 


(13) 


(14) 


(15) 


Here.  a(h)  is  the  fraction  of  ice  present  in  a  category  that  participates  in 
ridging,  k  is  a  constant  that  determines  the  ratio  of  final  to  original  thick 
ness  during  the  ridging  process  and  it  defines  the  fall  angle  of  the  ridge 
till  and  keel,  p '  is  the  coefficient  of  sliding  friction  of  ice,  and  if'  is 
th.  angle  of  friction  of  the  blocks  of  ice  In  the  ridge.  When  realist: 
constant  are  used,  friction  accounts  for  over  80"  of  the  energy  dissipat..’. 

W  .•>)  these  two  energy  sinks  are  combined,  the  large-scale  strength  becomes 


O' 


p*  *  p*  f  p* 

p  f 


p* 


a(h)  dh 


A 


1  ' 


g  k 


a 

t  an  <■ 


(16) 

(17) 

(18'' 


J 
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account,  the  noise  propagating  from  all  sources  to  any  site  may  then  be  simu¬ 
lated.  This  model  of  background  noise  is  based  on  the  fundamental  principles 
of  acoustic  propagation. 

Observed  background  noise  in  the  Beaufort  Sea  during  the  AIDJEX  field 
experiment  in  1975-76  is  used  to  test  the  model.  Low-frequency  noise  signals 
(10  and  32  Hz)  are  compared  with  the  simulated  noise.  Time  histories  of 
observed  and  simulated  noi6e  compare  well  in  frequency,  duration  and  magnitude. 
The  model  simulates  462  of  the  variance  of  observed  noise  during  a  120-day  test 
period  (the  correlation  coefficient  is  0.67).  The  comparison  is  broken  into 
20-day  intervals  to  determine  variations  in  time.  During  several  of  these 
20-day  intervals,  the  model  has  a  correlation  of  over  0.80,  explaining  over  642 
of  the  variance.  When  potential  error  sources  are  considered,  it  is  found  that 
much  more  of  the  variance  could  be  explained  by  the  model  if  several  error 
sources  could  be  reduced. 

Thfs  comparison  is  good  enough  to  warrant  further  study  and  testing  of  this 
new  concept.  First,  however,  data  sets  must  be  obtained  in  other  regions  and 
at  other  times,  since  various  noise-making  processes  will  contribute  different 
amounts  of  noise  in  different  areas,  such  as  at  an  ice  edge.  Also,  Ice  motion, 
wind  and  current  data  must  be  used  in  a  complete  tee  dynamics  simulation  model 
to  avoid  contamination  b>  velocity  inhomogeneities.  A  limited  amount  of  d > r a 
exists  for  areas  other  than  the  Beaufort  Sea,  and  recent  field  experiments, 

•uch  as  FRAM^  and  MIZEX,^  will  provide  more  data.  Further  testing  wit); 
these  additional  data  will  help  determine  the  range  of  appl ieabllity  of  th< 
concept  that  energy  dissipated  by  the  Ice  cover  Is  a  useful  measure  of  the 
sen  ..  .  *  jf  ambient  noise 
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The  statistical  approach  presented  here  Ignores  the  physical  causes  of  noise 
generation.  To  this  end,  this  paper  has  focused  on  showing  that  one  of  these, 
namely  the  power  dissipated  by  ridging,  is  highly  correlated  with  background 
noise.  Since  there  are  many  sources  of  noise,  we  cannot  expect  a  perfect 
correlation  with  only  one  source.  Instead,  the  simulation  model  can  be  used  to 
verify  qualitatively  the  contribution  of  each  noise-making  process  and  then  to 
determine  how  much  of  the  observed  noise  is  accounted  for  by  each  process. 
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Relative  intensity  of  acoustic  waves.  The  transmission  loss 
function  is  given  for  long-range,  low-frequency  propagation  over 
deep  water. 

Spatial  distribution  of  AIDJEX  sensors  in  the  Beaufort  Sea  on 
Julian  day  381  (January  16,  1976). 

Comparison  between  observed  background  and  simulated  noise  by 
ridging  for  120  days  in  winter  1975-76.  Data  are  presented  for  the 
one-third  octave  band  centered  around  10  Hz.  Sound  intensity  is 
presented.  Simulated  sound  intensity  16  nondlmensional  with  an 
arbitrary  scale  because  the  signature  function  S(f)  is  unknown. 

Comparison  between  observed  background  and  simulated  noise  by 
ridging  for  120  days  in  winter  1975-76.  Data  are  presented  for 
one-third  octave  bands  centered  around  32  Hz.  Sound  intensity  is 
presented.  Simulated  sound  intensity  is  nondlmensional  with  an 
arbitrary  scale  because  the  signature  function  S(f)  is  unknown. 

Observed  noise  and  simulated  noise  from  energy  dissipated  by 
ridging,  qr,  at  10  Hz, 

Observed  noise  and  simulated  noise  from  energy  dissipated  by 
ridging,  q^,  at  32  Hz, 
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TABLE  I.  Transmission  loss  law  constants. 


Path  Length 


Less  than  185  km 
(100  nautical  miles) 


Greater  than  185  km 
(100  nautical  miles) 


Constant 

a 

b 

r 


69.3 

0.07 

-0.00081 


63.2 

0.03 

0.035 


d 


0.00026 


0.00059 
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TABLE  II.  Correlation  between  observed  background  noise  Intensity  and  noise 
Intensity  simulated  by  power  dissipated  by  sea  Ice  cover  in  the  central 
Beaufort  Sea  during  the  AIDJEX  main  experiment  (winter  1975-76).  All  data 
were  filtered  to  remove  frequency  content  higher  than  1  cycle  per  day  and  then 
aampled  at  3-hr  Intervals. 


Date  Correlation  coefficient,  r 

(Julian  day) 


RELATIVE  INTENSITY  TRANSMITTED 


•  2  Spatial  Distr  ibut  ion  oi  AIlMhX  Sensor--  ip  tlir  Bt  nf  m(  V 
on  Julian  -i a v  181  (  l.inn  irv  16,  19761. 


SIMULATED  1 0  Hz  SIMULATED  1 0  Hz  OBSERVED 

INTENSITY  o.(  x  10-®  iNTENSITV,  q  x  lO*8  10  Hz  INTENSITY  x  10"9 


OBSERVED  10  Hz  INTENSITY  x  10"9 


SIMULATED  32  Hz  INTENSITY,  qr 


Obsi  rvc-il  N.  >  j  am!  Si  1 1  .• r .  ■  !  Nr .  •< 
Ridj’lm’ ,  cl  a  •  12  lb- 


